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We use a single site lattice in four dimensions to study the scaling of large TV Yang-Mills field 
coupled to a single massless Dirac fermion in the adjoint representation. We use the location of 
the strong to weak coupling transition defined through the eigenvalues of the folded Wilson loop 
operator to set a scale. We do not observe perturbative scaling in the region studied in this paper. 
Instead, we observe that the scale changes very slowly with the bare coupling. The lowest eigenvalue 
of the overlap Dirac operator is another scale that shows similar behavior as a function of the lattice 
f*") coupling. We speculate that this behavior is due to the beta function appoaching close to a zero. 
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I. INTRODUCTION 



Particle accelerators experiments provide strict bounds for the beyond standard model physics. For technicolor 
Oh it means that the coupling constant has to exhibit walking behavior. Otherwise the theory cannot simultaneously 
explain the mass pattern of standard model fermions and the suppression of the flavor changing neutral currents 
[TH3] ■ Hence, lattice studies of vector like gauge theories with appropriate choice of fermion matter with the aim 
of understanding the conformal window has recently attracted considerable attention(see [5] and references therein). 
The gauge group is chosen to be SU(N) and the number and representation of fermions is such that the theory is 
expected to be conformal or near conformal [6]. 

£ L6t 

cn i 



P " cfN 

define the inverse 't Hooft coupling on the lattice. Let 



(1) 



t = lna (2) 

define the logarithm of a lattice scale where a(b) could be the square root of the string tension measured on the lattice 
at the coupling b. The beta function of the lattice is defined as 

m-<f. (3) 



The perturbative beta function leads off as 



P(b) = -6 - y + • • • . (4) 



As is well known [7J, only the one and two loop coefficients, bo and b±, are universal and the higher order coefficients 
in a Taylor expansion of /3(b) in powers of 6 _1 depend on the choice of a(b). In fact, one can imagine choosing an 
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a(b) such that all higher order coefficients are zero. We will not have such control on the choice of a(b). In particular, 
there is no reason to expect the location of the zero of the beta function to be independent of the choice of a(b) - all 
we can expect is for the zero to remain stable if it is at a perturbatively weak coupling. 

The choice of fermionic matter can be motivated by the presence of a zero in the two-loop perturbative beta 
function. In order to maintain asymptotic freedom, all choices are such that bo > 0. The two loop beta function has 
a zero if b\ < 0. Some of the choices currently being investigated are: 

• SU(3) gauge group with twelve Dirac flavors of fermions in the fundamental representation [81410] - b\ is negative 
if we have nine or more Dirac flavors but the zero of the two loop beta function occurs at smaller coupling for 
larger flavors. 

• SU(2) gauge group with two Dirac flavors of fermions in the adjoint representation [TUfTl] - This is the only 
choice based upon b and b\ since bo < if we choose three or more Dirac flavors and b\ > if we choose one 
Dirac flavor. 

• SU(3) gauge group with two Dirac flavors in the two-index symmetric representation [T5HT7] - In this case bi > 
if there is only one Dirac flavor. One can also choose three Dirac flavors and still maintain asymptotic freedom. 

The case of SU(N) gauge theory coupled to / flavors of Dirac fermions in the adjoint representation is interesting 
for two reasons: 

• The first two coefficients of the beta function are 

l 11-4/ , 17-16/ 

& 1 = mo_4 ■ (5) 
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and are independent of N [33] . The three interesting choices for a theory with an infra-red fixed point are 
/ = |, 2, | based on the two-loop beta function. 

Numerical evidence along with continuum arguments [IB]-[2S] suggest that Eguchi-Kawai reduction holds in the 
large N limit as long as one uses periodic boundary conditions for fermions. This is expected to be the case for 
/ > \ and for non-zero quark masses [19] . 

We have the possibility to study theories with an infra-red fixed point that have only four SU(N) degrees of freedom 
provided we consider the N — > oo limit. For finite N, the massless fermionic operator is a finite dimensional operator 
that decouples into chiral sectors. The fermion determinant is positive in each chiral sector and we can define a theory 
for any real value / since we can write 

(detp) f = e /lndot f>. (6) 

If Ijl > / > yj, the two loop beta function has a zero and the theory has an infra-red fixed point. 

Our aim in this paper is to use overlap fermions [28]- [30] and study the / = 1 theory on a single site lattice. We do 
not expect the beta function to have a zero from the perturbative viewpoint. Even if it has a zero, we expect it to be 
at strong coupling. With this in mind we expect a computation of the running coupling to agree with the two-loop 
running. Contrary to this expectation, we will show that the coupling runs much faster than what is expected from 
perturbative running. 

The model on the single site lattice is defined in Sec. [IT] We will numerically study this model using the Hybrid 
Monte Carlo (HMC) algorithm with pseudofermions as described in Sec. |III[ It is numerically difficult to extract the 
string tension. On the other hand there is an observable based on the Wilson loop operator [3TJ [35] that shows a 
transition from weak to strong coupling and we will use the location of this transition to set our scale as discussed 
in Sec. |IV A[ We will also look at the eigenvalues closest to zero of the overlap Dirac operator. We will compare the 
behavior of the lowest positive eigenvalue as a function of the lattice coupling and compare its behavior to the scale 
set using the Wilson loop operator. 

Results for the behavior of the scales set using the Wilson loop operator and the lowest positive eigenvalue of 
the massless Dirac operator are discussed in detail for the case of theory with massless fermions in Sec. [V] We will 
show that both scales are monotonic in the coupling and that they both vary very slowly with the coupling. We will 



speculate on the possibility of a near-zero of the beta function in Sec. VI 
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II. THE MODEL 



The action on a single site lattice with one flavor of adjoint Dirac overlap fcrmion is given by 

S = S g + S f . (7) 

The gauge action is 

1 4 

S g = -12bNP; P=- £ ^ [*W#^] . (8) 

where the four gauge degrees of freedom, (fi = 1, 2, 3, 4), are SU(iV) matrices. The lattice gauge coupling constant 
is b = . The overlap fermion action is 

S f = -flndetH (fi), (9) 
where the Hermitian massive overlap Dirac operator is defined by 

= \ [(1 + AO 75 + (1 - M) e(ff)] , (10) 

with n € [0, 1] being the bare mass. We note that the eigenvalues of H o (0) are in the range [—1,1] with exact zero 
eigenvalues and exact ±1 eigenvalues corresponding to a gauge background with non-zero topology. The Hermitian 
Wilson Dirac operator for adjoint fermions is given by 



H = 



4 - m - \ E M (V, + V*) \ (V, - V*) 

-\Y,A (V, - V*) -4 + m + \ E„ (V, + V*) 

(4-m) 7 B-EK^ + «'M. ( n ) 



where 



--ite -'") <i2) 

and y M are the link matrices in adjoint representation. The action of on $ is given by 

V^ = U^Ul; Vfi = Ul$U ll . (13) 
One can verify that H is Hermitian in the usual sense: 

Tr&HQ = [Tr$ t i/*] * = Tr t$] . (14) 

Therefore = (i/*) t and in addition it is also true that TriJ$ = if Tr$ = 0. The same is also true for H {p). 



III. THE NUMERICAL ALGORITHM 



We used the Hybrid Monte Carlo (HMC) algorithm to generate according to the measure 

Z = J [dU^e- 3 . (15) 

Let us introduce a Hamiltonian 

H=\Y,TtHl + S, (16) 

where matrices, fj, = 1,2,3,4 are elements of the su(N) algebra and conjugate to U^. The HMC algorithm 
involves the computation of the force, The gauge part of the force is simple to compute numerically, but the 
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fermionic part of the force is computationally intensive. An exact algorithm was developed in [55] to compute the 
fermionic part of the force. This algorithm scales like TV 6 . In addition to using this algorithm, we also developed 
a pseudo-fermion algorithm in order to compute the fermionic part of the force which scales like TV 4 . We present 
the details of the pseudo-fermion algorithm in this section. Both algorithms were used to obtain the numerical data 
presented in this paper. 
We note that 



(17) 



and 

detH (n) = detir£f(/i) = detH 2 _(n), 
in the zero topological sector. |34j The overlap fermion action can be rewritten as 



Sf = Tr 



(18) 



(19) 



where the pseudofermions have positive chirality and are traceless N x N complex matrices with an additional 
two component spin index. 

For numerical purposes, we will represent e(H) as 



«*0 = E 



k=l 



< Pi < p 2 ■ ■ ■ < p n , 



(20) 



with n chosen such that the representation is accurate in the spectral range of H 2 assuming some lower bound on the 
spectrum of H 2 . 

The algorithm starts with one choice for U^. Then, we draw according to a Gaussian distribution. We also 
draw Dirac indexed traceless Hermitian matrices, 'J/ according to the Gaussian distribution, Tr^'I', and set 



The equations of motion for are 



Setting £p = results in 



$+ = P + H {jj)V. 



dr dr 



0, 



and 



dSg_ 
dr 



E Tr [ H » D i\ 



^ = E^[^] 



The equation of motion for is given by 



dH. 



dr " " 



Taking the derivative of 5 S in d8h with respect to r and using ( 22 ) we arrive at 



(21) 



(22) 



(23) 



(24) 



(25) 



The derivative of Sf in (19 1 with respect to r using (17) is 



dr 



1-p 



Tr 



de(H) . 
dr 



r + = [H 2 + ((,)] 



(26) 



(27) 
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Substituting the representation (20 1 for e(H), we can write 

dS f _ 1 - m 2 — 



dr 



k=l 



where 



Using (11), (13) and (22), we can show that 



^2 ( HcPfcTr 



1 



dH . 



H 2 + Pk 



7-fcTr 



A'- 



tiff 



" fe dr " h 



(28) 



(29) 



Tr 



X^X 
dr 



J2^-[HMX)]. 



(30) 



where 



4 



t y 



xj^XjUl 



(31) 



for any complex matrix X. It is clear that At(X) = A fJi (X) and that TrA fl (X) = 0. Therefore, 



^ [rkPkA^Tk) - rkA^k)] 



(32) 



fc=i 



Given <I> + in (21 1, we compute T + in (27) with the standard conjugate gradient algorithm. Each action of H% + (/j,) 
that is part of the conjugate gradient algorithm involves the action of e(H) on a Dirac indexed traceless Hermitian 
matrix. We use the multiple mass conjugate algorithm for each action of e(H) represented by (20). The core of the 



fermion algo rithm is the action of H on a Dirac indexed traceless Hermitian matrix and this operation scales like 
N 3 [see (13)]. In addition, the computational cost depends on the gap of H and H 0+ (fi). The former is large and 
therefore does not seriously affect the computational cost. Since we are interested in studying chiral properties of the 
theory and want to work with as small a bare mass, //, as possible the smallest eigenvalues of H o+ (0) will scale like 
7V~ 2 and the condition number grows like iV 2 . 



IV. OPERATORS 



We will focus on measuring two quantities that will help us understand the running of the coupling with the scale 
and reveal numerical evidence for a singular point. One observable looks at the property of the gauge field and the 
other looks at the property of the massless fermion. 



A. Weak to strong coupling transition | 32| 

A folded L x L square Wilson loop operator in the fi — v plane is given by 

W(L) = U^Ul L Ut L . (33) 

The eigenvalues, e l6k ; k = 1, •• •, N of this operator are gauge invariant. Let p(6;L,b) be the distribution of these 
eigenvalues with 8 € [— tt, 7t). This distribution undergoes a transition at N — > oo as the size, L, is changed at a fixed 
coupling b: the distribution has a gap at 7r for small areas and it becomes gapless for areas larger than a critical area 
A c (b). There is a universal function describing the distribution in terms of the scaled variables derived from A(b) and 
9 in the vicinity of A c (b) and n. 
Let 



N (£;L,b) = (|dct(e5 + e -*W(L))J>; £ = tanh|. 



(34) 
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The region close to £ = probes 9 close to 7r. Let 

N (t L, b) = C (L, b, N) + C\(L, b, N)e + C 2 (L, b, N)^ + ■■■. (35) 
It is useful to define a Binder cumulant type quantity 

One can show using the universal scaling function that 

n(L c (b),b, oo) = Ao y i J = 0.364739936 (37) 
487r 

We can define L c (b, N) at a fixed A" and b as the length where 

n(L c (b, N),b, N) = 0.364739936, (38) 

and 

lim L c {b,N) = L c {b), (39) 

N— s-oo 

will be the location of the transition at infinite N. 

Since we are working at a fixed but large ./V in this paper, we will define our length scale as 

< b) = L^NY (40) 

B. Low lying fermion eigenvalues 

The eigenvalues of the massless Hermitian overlap Dirac operator, H o (0), can be used to see how they scale and 
if they show evidence for chiral symmetry breaking. The eigenvalues come in doubly degenerate pairs and there is 
also a pairing of positive and negative eigenvalues due to exact chiral symmetry on the lattice. We computed all the 
eigenvalues of the massless overlap Dirac operator. 

Let 0<Afc<l,fe = l,---, A^ 2 — 1 with < A^+i denote all the positive distinct eigenvalues where each eigenvalue 
is doubly degenerate and each positive eigenvalue has a negative eigenvalue pair. We can use 

A(6) = (AO (41) 

as another choice for our length scale. 

If chiral symmetry is broken, the chiral condensate sets a scale. In particular, we expect 

to be independent of the coupling for a few low values of k. As N increases, we expect more to be independent of 
the coupling. In addition, we expect X(b) to approach a finite limit as N — > oo. 

V. SINGLE SITE MODEL WITH MASSLESS ADJOINT FERMIONS 

Our choice of b and N are based on numerical feasibility. We expect the approach to the large Af limit to get 
slower as we increase b. Since the numerical costs increase rapidly with N, we cannot make N as large as we wish. 
Computational costs are manageable if we choose N in the range of 13 to 25. We will restrict ourselves to mainly a 
single value of N, namely, A^ = 18 and also provide some additional data with A^ = 25 to understand finite A^ effects. 
We have chosen the couplings in the range b £ [0.32, 0.70]. Our definition of the coupling is related to the conventional 
lattice coupling by 

P = 2bN 2 . (43) 
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N b 
18 0.32 
18 0.35 
18 0.40 
18 0.45 
18 0.50 
18 0.53 
18 0.55 
18 0.57 
18 0.60 
18 0.65 
18 0.70 
25 0.40 
25 0.45 
25 0.50 
25 0.55 
25 0.60 
25 0.65 



(P) 

0.6092(7 
0.6290(7 
0.6720(7 
0.7045(6 
0.7325(5 
0.7468(5 
0.7562(5 
0.7650(5 
0.7775(5 
0.7943(5 
0.8076(4 
0.6888(5 
0.7187(5 
0.7444(5 
0.7678(5 
0.7850(4 
0.8024(4 



0.4442(24) 

0.4251(23) 

0.3858(22) 

0.3561(29) 

0.3354(27) 

0.3050(36) 

0.2931(29) 

0.2820(28) 

0.2704(26) 

0.2566(29) 

0.2354(131) 

0.3761(27) 

0.3400(34) 

0.3134(18) 

0.2973(17) 

0.2864(15) 

0.2778(17) 



NX(b) 

0.0544(8 

0.0507(7 

0.0440(6 

0.0381(5 

0.0330(4 

0.0301(4 

0.0279(4 

0.0278(4 

0.0256(4 

0.0229(3 

0.0213(3 

0.0404(6 

0.0362(6 

0.0318(5 

0.0277(4 

0.0238(4 

0.0218(4 



TABLE I: Table showing the various values of couplings where simulations were performed with massless fermions with N = 18 
and N = 25. The results for the average plaquette and the two different choices for the scales are also shown. 



Our range of coupling corresponds to f3 G [2.56, 5.6] for SU(2) and (3 G [5.76, 12.6] for SU(3). The choice of couplings 
falls in the range of recent simulations with adjoint fermions using SU(2) as the gauge group [TTHT3"] and also in 
simulations with fermions in the symmetric two-index representation and SU(3) as the gauge group |15) . 

In addition to these physical parameters, we also have to choose the value of the Wilson mass parameter, m, in 



(111. It is an irrelevant parameter but needs to be chosen in a specific range to realize the correct continuum limit. 
Based on previous studies on a single site model with adjoint fermions 119) . we set m — 4 in this paper. 

Table [I] shows the various values of couplings where simulations were performed along with the results for the 
average plaquette (c.f. (|8|), a(b) (c.f. (40 1), and NX(b) (c.f. 41 1). A plot of the average plaquette is shown in Fig. T] 
The plaquette leads off as 1 — + 0(fr~ 2 ) where the coefficient of b^ 1 is not affected by fermions. A fit of the 
data shows a smooth approach to unity as b — > oo. The data also shows a measure of the fact that the eigenvalues 
of Polyakov loop operators Up are uniformly distributed. The four data points shown by different colored squares, 
correspond to the average values of 

^ = ^(l-^|Trtg 2 ); (44) 

for = 1, • • • ,4 with Pi < P 2 < P% < P4 on every gauge field configuration. An average value of | in the large N 
limit shows that the Zn symmetries are not broken. Our results are very close to \ for both N = 18 and N = 25 and 
we can assume that reduction to a single site holds and we are simulating an infinite volume theory. 
We define 



'tad 



HP), (45) 



as the tadpole improved coupling and plot the running of this coupling versus our two logarithmic scales, lna(6) and 
hi(NX(b)) in Fig.[2]and Fig. [^respectively for the data points listed in Table]]] The data with errorbars are shown with 
solid circles in both figures. We chose one point in the middle of the range as our renormalization point and the solid 
curve represents the result based on two loop perturbation theory. Clearly, there is no agreement and the coupling 
runs much faster than what is expected from two loop perturbation theory. This indicates that we are working with 
lattice couplings that should be considered as strong in spite of the fact that we used values that would be considered 
as weak in theories that do not have additional fixed points. 

We end this section by presenting some details pertaining to the two scales. We plot £l(L, b, N) for N = 18 and 
N = 25 in Fig. |4] We have set the x-axis to L/L c (b) where L c (b) is obtained using (38). We see that the value for 



f2 flattens out for large loops and this is a finite N effect. In addition, it flattens out at a higher value for weaker 
coupling and this is because finite N effects sets in at smaller physical loop sizes at weaker coupling. Due to these 




FIG. 1: Average value of the plaquette along with the average values for the four different Polyakov loop observables for 
massless fermions at N = 18. 



two effects, the determination of L c (b) at weaker coupling has larger finite N effects. We note that the finite N effect 
get weaker at N = 25 where we can perform a better estimate of L c {b) even at b = 0.65. The larger finite N effect at 
weaker coupling can also be seen in Fig. [2] 

The complete spectrum of the distinct eigenvalues of the massless adjoint overlap Dirac operator are shown in Fig. [5] 
for three different couplings at N = 18 and N — 25. All three plots show the same qualitative behavior. We see a 
concentration of small eigenvalues (less than 0.1) followed by a bulk like distribution. We think the distribution for 
A < 0.1 is due to the would be zero modes in the gauge field background that is diagonal. We believe that this part 
of the distribution will shows signs of chiral symmetry breaking if one exists. If this is the case, we would expect the 
lowest eigenvalue to scale like -^j. But we only see evidence for scaling like in Fig. [3] where a plot with ln(7VA(&)) 
on the x-axis show good agreement between N = 18 and N — 25. Therefore, we do not yet see evidence for chiral 
symmetry breaking. 

The plot of rfc as defined in (42) versus k is shown in a log-log plot in Fig. [6] Here again, one sees a separation 
between the low eigenvalues (the would-be zero modes in a diagonal background) and the bulk. Furthermore, the 
ratios do not change much with coupling for k < 5 for N — 18 and k < 7 for N = 25 showing that the finite N effect 
is of order ^ . 
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-1.5 -1.4 -1.3 -1.2 -1.1 -1 -0.9 -0.8 

ln(a(b)) 



FIG. 2: Running of the tadpole improved coupling versus the logarithmic scale lna(6) for massless fermions at N = 18 and 
A = 25. 



VI. SPECULATIONS 



How does one see the effect of an infra-red fixed point in a lattice computation? According to the two loop beta 
function, the zero occurs at 

bi 1 16/ -17 
Kif) ^~V = 8^^4T- (46) 

If / is close to the upper limit of the zero occurs at a value of the coupling that can be considered as perturbative. 
Different choices for a(b) will not matter in a study of the infra-red fixed point. One should see numerical evidence 
for a very rapid change of a(b) versus b close to &*(/) indicative of a zero of the beta function. The critical values 
for / = |,2 and | are 6* (|) = 6*(2) = ^ respectively. 6* (|) = None of these would be considered 

perturbative and it is quite likely that the study of the infra-red fixed point is strongly affected by the choice of a(b). 
In particular, it is quite possible that the location of the infra-red fixed point depends on the choice of a(b) and there 
is even a possibility that the existence of an infra-red fixed point on either side of the lower boundary of / = j| 
depends on the choice of a(b). We have set / = 1 in this study and we clearly do not see perturbative behavior as 
discussed in Sec. W\ 
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FIG. 3: Running of the tadpole improved coupling versus the logarithmic scale lnA(&) for massless fermions at N = 18 and 
AT = 25. 



In order to understand the behavior of the running coupling better, we consider beta functions of the form 



A: 13(b) 

B : p{b) 

C : /3(b) 

D : 13(b) 



-b (1 - 
-b (l - 



1 - 



-b 



1 



-&o 1 1 - 



-Mii-Vr 



a 

for b > 6* 
for b < b* 



(47) 



where it is assumed that all parameters except s are positive and we also assume that bj > bjj. The parameter s 
could be positive or negative. 



• Case A is the two loop beta function written for the case where it has a zero, bj is the location of the infra-red 
fixed point. The beta function is shown in solid red in Fig. [7] and the a plot of the coupling with the scale is 
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FIG. 4: Plot of the quantity, Q(L, N,b), as a function of L/L c (b) at iV = 18 and iV = 25 for massless fermions at several 
different coupling. 
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FIG. 5: The full distribution of the eigenvalues of the massless adjoint overlap Dirac operator for three different couplings at 
N = 18 and N = 25. 



shown in solid red in Fig. [8j The logarithmic scale goes to positive infinity as b — > bj . In order to define the 
continuum limit, we need to take the limit t — > —oo and this is achieved by taking b — > oo in the usual manner. 
This is what one expects to see if nothing occurs beyond what is expected in two-loop perturbation theory. 

• Cases B and C have an additional p- term and were discussed in |27j . 

— The beta function and the dependence of coupling on the scale for case B are shown in dashed green in 
Fig. [7] and Fig. [8] respectively. The logarithmic scale goes to positive infinity at b = bj and goes to negative 
infinity at b = bjj which is an ultra-violet fixed point. In addition to defining a continuum limit as b — > oo, 
we can also define a continuum theory as b — > bjj. It is likely that the location of b\j depends on the choice 
of a(b). If we assume that a bjj exists for all choices of a(b), then we can define a continuum theory by 
taking a limit &—>&[/ that is operator dependent. 

— The case B with bj = bjj = 6* is same as case C with a = 0. The beta function and the dependence of 
coupling on the scale for this case are shown in dot-dot-dashed brown in Fig. [7] and Fig. [8] respectively. 
The limit of b — > 6* from above results in the logarithmic scale approaching positive infinity and the limit 
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FIG. 6: The ratios of the eigenvalues of the massless adjoint overlap Dirac operator are shown in log-log plot for all couplings 
common to TV" = 18 and N = 25 in Table. U 



of 6 — > 6* from below results in the logarithmic scale approaching negative infinity. A continuum theory 
can be defined in the limit of 6 — > 6* from below. As before, 6* is expected to depend on the choice of a(b). 

— The beta function and the dependence of coupling on the scale for case C with a > are shown in dot- 
dashed blue in Fig. [7] and Fig. [8] respectively. In this case, there is a region around 6 = 6* where the 
coupling walks. How far it walks, depends on a and the range [ti,tu], where the coupling essentially 
remains a constant, grows as a decreases. As before, 6* is expected to depend on the choice of a(b). The 
continuum limit is defined as 6 — > oo. 

• In all cases discussed above, the location of 6/ and/or bjj depends on the choice of a(6). Case D is different 
in this respect. If the beta function has a zero that is non-analytic, we expect the location of the zero to not 
depend on a(6). In other words, if one observable shows non-analytic behavior at some coupling, all observables 
are expected to show non-analyticity at the same coupling. If s = —1 and p + = p_ = 1, we recover case A. If 
s = 1 and p+ = p— =2, we recover case C with a = 0. The beta function and the dependence of coupling on 
the scale for s = 1 and p+ = p- = | are shown in dot-dash-dashed orange in Fig.[7[and Fig. [8] respectively. Like 
in case C, for this choice of parameters, there is a region around 6 = 6* where the coupling walks. Whether 
there is a infra-red/ultra- violet fixed point at 6 = 6* depends on the choice of p± as can be seen from Fig. [9] 

— If p- < 1, the scale change from 6 = to 6 = 6* is finite. In such a situation, we cannot define a continuum 
theory by taking 6 — > 6* from below. 

— If p- > 1, the logarithmic scale goes to negative infinity as 6 — > 6* from below and 6* is an ultra-violet 
fixed point from below. We can define a continuum theory in this limit. 

— if p + < 1, the logarithmic scale starts out at a finite value for 6 = 6* and goes down to negative infinity as 
6 — > oo. A continuum theory can be defined as 6 — > oo. 

— If p+ > 1, 6 = 6» is aji infra-red fixed point from above and the logarithmic scale goes to positive infinity 
as 6 — » 6* from above. A continuum theory can be defined as 6 — > oo. 

As the maximum of /3(b) moves away from zero in case C, the range where the coupling runs slowly shrinks. As 
the power, p±, goes below unity in case D. the range where the coupling runs slowly shrinks. In such cases, it will 
be difficult to see evidence for a slow running of the coupling. Instead, one will see a fast running of the coupling 
on either side of the maximum in case C or the location of the singular point in case D. The numerical data shown 
in Fig. [2] and Fig. [3] is close to what is seen in case C or case D for couplings around 6* . Since we see do not see 
agreement with perturbation theory and we see a running that is significantly faster, the couplings we are using on 
the lattice are probably on either side of 6* and / = 1 probably corresponds to a case where the region of slow running 
of coupling is very small. Clearly, our data is not precise enough to distinguish between either of these cases. But our 
data suggests the presence of a 6* . 
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FIG. 7: Different cases for the beta function are shown. 



VII. CONCLUSIONS 



The single site model of a large N gauge theory coupled to massless adjoint fermions was numerically studied in 
this paper. We have studied this model with a single flavor of adjoint fermion numerically using the Hybrid Monte 
Carlo algorithm and pseudofermions. We studied the running coupling using two different choice of scales and did not 
find agreement with two-loop perturbation theory at intermediate values of the tadpole improved coupling. The two 
different choices of scales were the transition from weak to strong coupling and the lowest eigenvalue of the massless 
overlap Dirac operator. This is the main result of our paper. 

Since one flavor of adjoint fermion is close to j| which is the lower bound for the number of flavors for a perturbative 
zero of the beta function, we speculate on the possibility that a near zero of the beta function might be the cause for 
our result being not in agreement with perturbation theory. Our numerical data cannot show with definiteness that 
there is a near-zero of the beta function but the behavior suggests such a possibility. 

The results presented in this paper are exploratory in nature and future simulations with different values of / will 
give a clearer physics picture. However, the work lays the foundation for the careful study of ultra-violet /infra-red 
fixed points in matrix models that mimic large N gauge theories coupled to adjoint fermions. We have the ability 
to treat the number of fermion flavors as a real number in the matrix model and study the presence of singular 
behavior in the associated beta function. It is likely that the behavior at the singular point, results in it being an 
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FIG. 8: The running of the coupling with the scale is shown for the different cases of the beta function. 



ultra- violet/infra- red fixed point for some range of fermion flavors. The numerical procedure developed in this paper 
for the case of a single Dirac flavor paves the way for future numerical studies of the matrix model with varying 
number of flavors. 



Acknowledgments 

R.N. acknowledges partial support by the NSF under grant number PHY-0854744. R.N. would like to thank Erich 
Poppitz for several useful discussions. The numerical calculations presented in this work have been performed on 
the Horseshoe6 cluster at the University of Southern Denmark (SDU) funded by the Danish Center for Scientific 
Computing for the project "Origin of Mass" 2009/2010. 



[1] B. Holdom, Phys. Rev. D 24, 1441 (1981). 

[2] K. Yamawaki, M. Bando and K. -i. Matumoto, Phys. Rev. Lett. 56, 1335 (1986). 



15 




t 



FIG. 9: The running of the coupling with the scale is shown for case D with different choices for p± with s = 1 and bo = 1. 



T. W. Appelquist, D. Karabali and L. C. R. Wijewardhana, Phys. Rev. Lett. 57, 957 (1986). 

J. R. Andersen, O. Antipin, G. Azuelos, L. Del Debbio, E. Del Nobile, S. Di Chiara, T. Hapola, M. Jarvinen et al, 
[arXiv: 1104.12551 [hep-ph]] . 

L. Del Debbio, PoS LATTICE2010, 004 (2010). 

F. Sannino and K. Tuominen, Phys. Rev. D 71, 051901 (2005) hep-ph/0405209 . D. D. Dietrich, F. Sannino, Phys. Rev. 
D75, 085018 (2007). |hep-ph/06 11341] . 
S. Weinberg, C ambridge, UK: U niv. Pr. (1996) 489 p 
A. Hasenfratz, [arXiv:1106.5293| [hep-lat]]. 

T. Appelquist, G. T. Fleming, M. Lin, E. T. Neil, D. A. Schaich, [arXiv:1106.l2T48| [hep-lat]]. 



Z. Fodor, K. Holland, J. Kuti, D. Nogradi, C. Schroeder, Phys. Lett. B703, 348-358 (2011). arXiv:1104.3124 [hep-lat]]. 

S. Catterall, L. Del Debbio, J. Giedt, L. Keegan, |arXiv:1108.3794| [hep- phi], 

T. DeGrand, Y. Shamir, B. Svetitsky, Phys. Rev. D83, 074507 (2011). [arXiv: 1102 1.2843] [hep-lat]]. 

A. J. Hietanen, K. Rummukainen, K. Tuominen, Phys. Rev. D80, 094504 (2009) .~ |a7Xiv:0904.0864| [hep-lat]] . 

F. Bursa, L. Del Debbio, D. Henty, E. Kerrane, B. Lucini, A. Patella, C. Pica and T. Pickup et al, Phys. Rev. D 84, 

034506 (2011) | arXiv: 1104.4301] [hep-lat]] . 

Y. Shamir, B. Svetitsky, T. DeGrand, Phys. Rev. D78, 031502 (20 08). |arXiv:0803 1707] [hep-lat]]. 
Z. Fodor, K. Holland, J. Kut i, D. Nogradi, C. Schroeder, [arXiv:1103.5998| [hep-latjj. 
J. B. Kogut, D. K. Sinclair, arXiv:1105.3749l [hep-lat]]. 



16 



Bringoltz, M. Koren, S. R. Sharpe, [arXiv:1106.5538| [hep-lat]]. 

Hietanen, R. Narayanan, Phys. Lett. B698, 171-174 (2011). arXiv:1011.2150| [hep-lat]]. 



B. 

A 

S. Catterall, R. Galvez, M. Unsal, JHEP 1008, 010 (2010). [arXiv: 1 006 . 2469| [hep- lat] ] . 

T. Azeyanagi, M. Hanada, M. Unsal, R. Yacoby, Phys . Rev. D82, 1250 13~f2010). [arXiv:1006.0717| [hep-th]] . 

A. Hietanen, R. Narayanan, JHEP 1001, 079 ( 2010). [arXiv:0911.2449| [hep-lat]]. 
E. Poppitz, M. Unsal, JHEP 1001, 09 8 (2010). |arXiv:091 1.0358 1 [hep-th]]. 

B. Bringoltz, JHEP 0906, 091 (2009). [arXi v:0905.2406 [ hep-lat |[. 

G. Cossu and M. D'Elia, JHEP 0907, 048 (2009) arXiv:0904.1353 [hep-lat]]. 

P. Kovtun, M. Unsal, L. G. Yaffe, JHEP 0706, 019 (2007). |hep-th/070202l| [HEP-TH]]. 

D. B. Kaplan, J. -W. Lee, D. T. Son, M. A. Stephanov, Phys. Rev. D80, 125005 (2009 ). |arXiv:0905.4752] [hep-th]]. 
R. G. Edwards, U. M. Heller, R. Narayanan, Phys. Rev. D59, 094510 (1999). |hep-lat/9811030| . 

H. Neuberger, Phys. Lett. B417, 141-144 (1998). [hep-lat /9707022]. 

R. Narayanan, H. Neuberger, Nucl. Phys. B443, 305-3 85 (1995). |hep-th/9411108| . 
R. Narayanan, H. Neuberger, JHEP 0603, 064 (2006). |hep-th/ 0601210 . 

R. Narayanan, H. Neuberger, JHEP 0712, 066 (2007). |arXiv:0711.455l| [hep-th]]. 

Trivial coefficients of N get absorbed since we have used the 't Hooft coupling and not g 2 . 

We are assuming that global topology is completely suppressed and one can restrict the theory to the zero topological 
sector. 



